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Abstract. We compute the production of particles from the gravitational field of 
an expanding mass shell. Contrary to the situation of Hawking radiation and the 
' production of cosmological perturbations during cosmological inflation, the example of 

an expanding mass shell has no horizon and no singularity. We apply the method of 
'ray-tracing', first introduced by Hawking, and calculate the energy spectrum of the 
produced particles. The result depends on three parameters: the expansion velocity 
(N ', of the mass shell, its radius, and its mass. Contrary to the situation of a collapsing 

mass shell, the energy spectrum is non-thermal. Invoking time reversal we reproduce 



o 

Hawking's thermal spectrum in a certain limit. 

a -1 

bfj. 

1. Introduction 

X 

Gradients in a gravitational field may cause the production of particles out of the 
vacuum; this is similar to the situation in quantum electrodynamics (QED). In the latter 
case a virtual pair of particles is ripped apart by the electric field. In contrast to QED, 
the gravitational 'charge' — the mass — is always positive; in this case, the tidal forces 
are thus responsible for the separation of the particle pair. But there is also another 
difference. The energy that is needed for the production of particles in a charged vacuum 
is taken out of the field energy. Since the energy of the gravitational field cannot be 
localised, the energy gain must be explained otherwise. One possibility is the existence 
of a horizon. Another possibility is the time dependence of the gravitational field. 

We compute the particle production in a time-dependent gravitational field induced 
by an expanding mass shell. This is done at a semiclassical level. The gravitational field, 
which acts through curvature of space-time, is treated as a classical background field in 
which a quantum field propagates. We consider a massless scalar field. So it is quantum 
field theory in curved space that we need to apply. 

Perhaps the best known result in this field is the thermal radiation emitted by a 
black hole. This effect, discovered by Hawking Q, shows a deep connection between 
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gravity, thermodynamics and quantum physics. Hawking's result was reproduced in 
several ways, but also other problems such as Unruh's accelerated observer [[J or the 
gravitational analogue to the Casimir effect || were studied. 

In cosmological models, the production of particles plays a central role in explaining 
the origin of the structures observed in the Universe. During an epoch of cosmological 
inflation, quantum fluctuations of space-time and matter are amplified by the accelerated 
expansion of the Universe to become the seeds of large-scale structures Q . 

All mentioned examples involve either a horizon, a boundary, or a singularity. To 
our knowledge, there exists no example in which particle production in an infinite space- 
time without horizon and without singularity has been calculated. The situation of an 
expanding mass shell represents such an example. We do not expect this work to have 
a direct astrophysical application, but the situation under study might be of interest 
when black holes evaporate. One could imagine that Hawking radiation is emitted 
from an evaporating black hole in a final burst, which we could then describe as a thin 
shell, carrying a mass somewhat larger than the Planck mass. If the emitted radiation 
consists of massive particles, this mass shell moves at a speed slightly below the speed 
of light and our calculation might be applicable. Since the expanding mass shell itself 
leads to gravitational particle production, this work might be useful for an estimate of 
backreaction effects in the final stages of black hole evaporation. We will not study this 
aspect here, the purpose of this paper being just to calculate the energy spectrum of 
the particles produced by an expanding mass shell. 

There are three main procedures that are used to deal with quantum fields in 
curved space-time: Hawking's classical ray-tracing JIJ, renormalization of the energy- 
momentum tensor (e.g. by point-splitting) || [J, particle production in the reference 
frame of a uniformly accelerated observer [EJ, [7| , from which we will choose the first one. 
A useful overview on computational matters can be found in || || . 

The procedure of ray-tracing in a semiclassical limit is appropriate if the energy 
carried by the quantum field is small with respect to the energy of the source of the 
gravitational field. Furthermore we have to stay on scales on which possible effects of 
quantum gravity do not occur. As Wheeler showed in the characteristic length for 

this to happen is the Planck length lp = \ ~ 10~ 33 cm. 

So we have to restrict our results to cases in which the effects of interest take place 
on scales ^> lp and ^> tp. However, it seems that there is no reason to restrict the 
calculation to masses M > mp, as long as the mass shell starts to expand at a radius 
3> lp. We set G = c = h = 1 in the following. 

This paper is organized as follows. First we will need to find an observable of particle 
density in curved space, which will be done in Section 2. After this we will introduce 
our model of the expanding mass shell and discuss its properties. The computation 
of the number density operator is presented in Section 4. The results are discussed in 
section 5, which contains plots of the spectral energy flux and total energy flux through 
a 2-dimensional surface. We conclude with a short summary. 
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2. Quantum fields in curved space-time 

Let us first recall some of the basics of quantum field theory in curved space-time 
||, H Physics in flat space-time is invariant under Poincare transformations. This 
singles out a set of (Fourier) modes and a vacuum state, which are invariant under 
Poincare transformations. In curved space-time there is no global Poincare symmetry. 

In the following we consider a minimally coupled scalar field ip. Its Fourier modes 
in flat space are attached to the energy uj by 

id t Uk = ujUk with uj > . (1) 

This uniqueness is lost when the coordiates are chosen arbitrarily. In curved space-time 
we have no distinguished time coordinate, thus no distinguished energy belonging to 
it, and no distinguished set of modes that can be related to the notion of particles. In 
general we have to deal with an arbitrary set of modes (ui) and ask for its relation to a 
different set (t>j). The quantum field <p can then be expanded in both sets: 

<p = ^(uidi + u*aj) , (p = ^{vjbj + v*bj) . (2) 

» 3 

The relations between these sets are given by the so-called Bogoliubov transformations 
[ TT| and have the form: 

V J = ^2( a Ji u i + PjiO , u i = ^2( a *ji v 3 ~ hrf) ( 3 ) 



with 

a jk = (vj,u k ) , Pjk= -(vj,u* k ) , (4) 
and the scalar product given by 

-i / [ Ui {d^) - (d^v*] v^i" d£" . (5) 

Here S is a space-like hypersurface with volume element dE, dll^ = nME, and n^ 1 is a 
normal vector of E with norm n^n^ = 1. 

The creation and annihilation operators (0) depend on the modes, and so does the 
vacuum state defined with respect to them. The vacuum state |0&) belongs to bi and 
is defined by [ 0&) = 0. In general, | Ob) is not identical with |0 o ), for which we have 
di\0 a ) = 0. It follows that the vacuum expectation value of the number density operator 
d\di must not necessarily vanish when applied to the vacuum |0&), 

{%\a\a i \%) = Ifyil 2 ■ (6) 

3 

This operator has a well defined meaning as the number (mode) density operator only 
if the modes used for the expansion can be put in a suitable relation to the known 
modes of Minkowski space-time to render an interpretation possible. This is possible if 
we are interested in the change of vacuum between two asymptotically flat regions of 
space-time. 
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If the eigenvalues are not discrete (z — > k), we have to deal with an integral that 
may have divergences. These divergences caused by the limit from periodic boundary 
conditions with discrete spectrum to an infinite volume V = (0 1 0) can be cured by 
normalization. This may be accomplished by dividing through the same divergence: 

(0 6 |a^j0 6 } ren = i(0 b |%|0 6 > = i y fe| 2 d 3 k' . (7) 

If (0fe|ata^|0fe) has the form 5(k — k)f(k, k), one finds 

(0 b \aLa^\0 b ) ren = — lim(0 b |atag|0 b ) = f(k) , (8) 
v k^k 

which will be useful later on. Now (0b|ata^|0ft) ren is the well defined spectral number 
(mode) density. 



3. The expanding mass shell 

We consider a spherical shell of mass M and radius R. Outside of the shell, space-time 
is described by the Schwarzschild metric; inside it space-time is flat. We assume that 
the thickness of the shell in negligible with respect to its radius (thin shell). Originally 
the shell's radius is constant, R = Rq- The shell is 'ignited' at some time r = (as 
measured in the rest frame of the total mass inside the shell). Afterwards the shell 
expands as R = R(t). The situation is illustrated in Fig. [j]. 

Our first task is to write down the metric of the space-time describing an expanding 
mass shell. The surface of the shell is the same whether it is measured from the inside or 
from the outside. We can thus choose the same radial coordinate r on both sides of the 
shell and we obviously can choose the same angular coordinates 6 and 0. Outside the 
shell we denote the time coordinate by t, which is the time measured by an asymptotic 
observer, who is at rest with respect to the centre of mass of the shell. The relation 
between the time coordinates t and r remains to be fixed by the boundary conditions on 
the mass shell. We take null coordinates and ignore the r 2 dfl 2 part of the line element 
for the moment. Inside the shell these coordinates are denoted by (U, V), and by (u, v) 
outside, where the line element is given by 

ds 2 = j(r) du dv , j(r) = 1 , (9) 

with 

f dr' f r dr' 

u + u :=t- —tt;' v + v :=t+ —r^- (10) 

Jr i\n J Ro 7(n 

Inside we simply have 

ds 2 = dUdV, (11) 

with 

V + V := r + r - R , U + U := r - r + R . (12) 
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Figure 1. The expanding mass shell (thick line). The thin lines represent the rays at 
the border lines between the three cases of ray-tracing. 



The outgoing light rays run on world lines with constant u and U. We fix the origin of 
the null coordinates by the convention that the incoming ray v = V = goes through 
the point (r = 0, Rq) and that U = u = denotes the very same ray on its way out (see 
Fig. [I]). Thus v = Vq = and at r = we have U + U = V + 2R . With our choice 
for U we finally have Uq = 2Rq. 

We define the expansion velocity of the shell by 

V '= df ' (13) 
which is actually the coordinate velocity as seen from the inside of the shell. In principle 
v is a function of time, but we will see below that it is a good approximation to work 
with a constant expansion velocity. In that case the world line of a point on the shell is 
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given by 



Ro + t^V = R + 1 ^(U + U ) when r>0 



R when r < 

Our result will finally depend on three parameters: Rq, M and v. 



(14) 



From the coordinate transformations (|T0|) and (|T2|) and from the definition of the 
expansion velocity (|i~3l) we find 



dw 

d!7 

df 
dV 



at 



r=R 



At j(R) )j i 1 ~ V ) Whel1 U > ~ U 

l/\/7o when U < —U 



dt i v 



where 



7o := 7(i? ) = 1 



dr 1 7(7?) 

1/V7b 

2M 



1 + v) when V > 
when V < 



R(\ 



(15) 
(16) 

(17) 



The matching of the inside and outside space-time demands the continuity of the line 
element 



(7 du dv)\ r=R 

which leads, with (|T5|) and 
dt 1 

d7 ~ ~f(R) 



(dU dV)\ r=R 
to 



7 (i2)(l 



1/" 



+ z/ 2 , 



(19) 



which describes the relation between t and r, and finally fixes the space-time geometry. 
Note that Eq. ([T^) describes nothing but the gravitational redshift z for a photon that 
is emitted inside the shell and crosses the shell at time r, i.e. dt/dr = 1 + 2. 

To trace a radial light ray through space-time, we need to find the transition 
functions between the two coordinate patches on the surface of the mass shell. The 
setting is illustrated in Fig. [I]. There are three different cases for a ray that arrives at 
infinity. Let us start with the trivial case (3). Everything happens before the onset of 
expansion. We do not expect any effect then. For the second case (2) the rays enter the 
shell before the start of the expansion and the shell is 'ignited' while the ray is crossing 
the interior. Only a final intervall in u of length uq is involved. For case (1), entrance 
and exit of the ray occur while the expansion is taking place. 

We obtain a collapsing shell by inverting the setting (u — ► —u, v — ► —v, r — > — r, 
etc.). To reproduce a collapse to a black hole, we need to take R — ► 2M and v — > lfj 
This procedure will allow us to reproduce Hawking's result from the setting of the 
expanding shell. 

The transition functions, which follow from (|I~5"D and ( |i~6"D by integration, are too 
complicated for an exact analytic treatment of the problem. We therefore have to find 

I After this change of the time direction the in- vacuum is defined in a different region, but the vacuum 
expectation value of the number density operator is still given by the same formula. 
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a reasonable approximation. Let us start from the inverse of Eqs. ( |T5| ) and ([H]), which 

reads 



dU 
du 
dV 
dv 



U>-U 



[7(1 - v 2 ) + u 2 }^ + v 

(1 + f) i 
[7(1 - z/ 2 ) + V 2 )? - V 



(20) 

(1-,) " (21) 

As a first approximation we assume that the expansion velocity v is constant. Secondly, 
we expect the biggest effect when the space-time curvature is high, thus in the vicinity 
of the 'explosion'. Consequently we expand the transition functions in the Newtonian 
potential difference — 2M(1/ R— 1/Rq), which is a small quantity at the beginning of the 
expansion. Case (3) is trivial. For case (2) we approximate dU/du around the 'ignition', 
R = Rq. For case (1) we have to consider two points at which the ray crosses the shell. 
The effect from the entrance of the ray into the shell will be larger than the effect from 
the exit of the ray. Thus we need to approximate dV/dv around 70 and dU/du around 
j(R(U = 0)) =: 7i. Around 70 we include terms up to linear order, whereas around 71 
we keep only the leading term. These expansions yield 

, 2M 

c — c n/rrN around 70 



where 



dU 




du 


u>-u { 


dU 




du 


[/>-£/( 


dV 




dv 


V>0 




dU 


c :-- 


du 




dU 


Ci :-- 


du 


do ■-- 


dV 


df 



70 



>R{U) 
i Ci around 71 

- - rfo^y 
+ c' (l -70) , 



around 70 



71 



in 



+ 4(1 -To) 



d_ dU 
c?7 du 



^dK 

d ° d 1 dv 



70 



7() 



As a consequence 
du 
dU 
dv 
dV 



Ro(l + u(U + U 



c [R (l -u) + u(U + U Q )] - c 2M(l - v) 
Ro(l + u) + uV 



(22) 
(23) 
(24) 

(25) 
(26) 
(27) 

(28) 
(29) 



d [Ro(l + v)+vV]-d' 2M(l + v) 
We denote the approximate transition functions around 70 at the entrance with 
£oOOj aroun d 70 at the exit with f] (u), and around 71 at the exit with r)i(u). The inverse 
functions of the last two are written as ^^(U) and ?7 ] ~ 1 (?7), respectively. Integration of 
(E3), (BH) and fl2D|) provides the functions that will be relevant to our calculation 



Vo 1 (U) = - + 

Co 

rAu) = -, 



U c 2M(l 



V) 



■In 



1 + 



CqUU 



(c R - c 2M)(l - v) + cquU _ 



,(30) 
(31) 
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«v, + 



(d R -d' 2M)(l 



(32) 



4. Particle production 



To compute the particle production caused by the moving mass shell, we consider 
an observer who defines vacuum before the onset of expansion (at J?~) and ask how 
this vacuum is seen by an observer at J^ + , i.e. we have to calculate the Bogoliubov 
transformations between the vacua at J?~ and J^ + . We can do this by following the 
track of the light-rays. 

We restrict our attention to a massless, minimally coupled scalar field, which obeys 
the equation of motion D(p = 0. We expect a damping of the amplitude with r~ l because 
of the spherical symmetry. Furthermore we split off an angular component in the form 
of spherical harmonics Yi m (9, 0): 

<p(r, t, 6,0) = J2 l ~Y lm (0, <t>)%{t, r) . (33) 

l,m 

With this ansatz, the wave equation in null coordinates becomes 

d 2 ^ l , , (2M l(l + l)\ T 

" l(r) -r + % ■ (34) 



dudv \ r 3 r 2 

The right-hand side has a minimum outside the Schwarzschild radius and acts like a 
potential well. Because of this, incoming rays may be reflected at the gravitational 
potential. To a good approximation one may neglect this effect of backscattering by 
simply dropping the right-hand side of (0). In the same way we neglect potential 
terms in the inside. The index I is left out in the following, since \I/ no longer depends 
on I if the centrifugal potential is neglected. 

We are searching for solutions of the equations 

outside, = inside. (35) 



dudv ' dU dV 

Any function that does not mix the two variables u, v or U, V respectively is a solution. 
We are looking for an incoming mode which turns into an outgoing one in a smooth 
manner and passes through r = 0. These solutions must behave like spherically 
symmetric Minkowski modes at J?~ . At r — ► oo the phase of an incoming wave is 
given by 

e~ iwv . (36) 

Inside of the mass shell, the corresponding mode travels on a curve with constant V. 
At r = we have U = V where the incoming ray exp(— iuV) turns smoothly into the 
outgoing ray exp(— iuU). At outgoing modes are proportional to 

e- iu ' u . (37) 

With the help of the transition functions we can now glue curves of constant v to 
curves of constant V and curves of constant U to curves of constant u. Figure |2] gives 
an overview of this procedure. Here we already inserted U$ = 2R Q . 
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v > 

-iwv >- g— iu}£(V) 



-iw£(»7(u)) 



case (1) 



« < o 



-iwV/^/To 



C/ > -2i? 

e -iuU/^o . e -iw V (u)/^o case (2) 



U < -2R 



-iw(u+«o-2/?o/ v ^o) case (3) 



Figure 2. Tracing of the incoming ray. 



In the following we restrict our analysis to the most interesting case (1), which 

provides the result for the extended expansion phase of the mass shell. Case (2) is 

restricted to a short time after the ignition of the shell, and thus is only characteristic of 

the ignition phase of our problem. Its importance depends on the velocity of the shell. 

So we now calculate the coefficients We expand the outgoing modes into plane 

waves. The outgoing modes are functions that depend on u only, and we obtain 

1 f°° 1 

: Al(u') =e- iw '" do/ for u > , (38) 



,-iutoim («)) 



47rV2c</ 



with 



a£l., u/ > 



(39) 



^ < 

This is a Fourier transformation, which allows us to compute the Bogoliubov coefficients 
without evaluatiing the scalar product ([5]). A Fourier transformation of the outgoing 
modes gives 

_ 1 



/ := 



■in 



-iw£o(»7i(«)) giw'w ^ u 



(40) 



Using the approximations for £o(^0 and rji(u), and introducing the notation A 
d' 2M(l + v)/(d%v), we find 



I = — 

2tt 



— i(£JjL— w')u 

e d o ' 



(Ro - d' 2M/d )(l + u)/v + ci« 



-iu;A 



(i2o-d(,2M/do)(l + i/)/i/ 
We further define x — x := Ciu/d and x := (i?o — d' 2M/do)(l + u)/(ud 
transition coefficients then become 



du . (41) 
The 



2n V oj 



o 



do 

Cl 



.r 



(42) 



where o;+ := a; + a;'—. 



The vacuum expectation value of the number density operator is given by 

2 



mj 2 du/ 



(t 1 







47T 2 cf 







0> 



da;' . 



(43) 
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This quantity is divergent and has to be regularized. First we will shift the u/-integration 
into the complex plane for some e. We do this because, for Imw + < 0, the inner 



integrations are well defined and may be expressed as incomplete T-functions [12 
Further we regularize the number density with the help of a 'frequency split' 

|/£,J 2 = nm/£,J?L* • (44) 

LO — >LJ 

Finally, we take the difference to the limit Xq ^ 0, which allows us to subtract the 
divergent part in a well defined manner. The point is that ([|3|) reduces in this limit to 
the expression for the number density that shows up in the Hawking effect 0, 

5{uj — uj) 



lim n / I Aid 2 du/ = lim 

xq^O 



cd^uj e 2?TwX — 1 

The only difference is that A replaces the inverse of the surface gravity, 1/k = 4M, that 
would show up in the situation of a collapsing mass shell. 

As a first step we rotate the integration variable x in the complex plane and obtain 



(45) 



|/tJ 2 d^ = hm 



o An 2 c\ 



oo-ie j 



0-ie 



e -iu + x x -iw\ dx 



du' . (46) 



Introducing uj + = uj + uj'— and uj := uj + e, we regularize the product of the inner 



integrations by a frequency split, 



f°° 1 
Jo IZ° 47r c l 



oo— le , ,/ i*xq~ 100 
UJ 



0-ie 



UJ 



-^+x x -^X dx 



J-t) 



xo+ioo 



(47) 



■t'o 



where we took the limit partially. In the first of the inner integrals we further substitute 
iuj + x =: y and in the second — i£u*,x =: y, 



|/3iJ 2 da/ = lim-^e-™ A > 



E ^o 4ti 2 c 2 uj 



(4? 



uj 

0-ie W + 



e -y y -iu\ dy 



Since 



ojl + e) 



K) 



-id; A 



0(e) 



e~ y y lwX dy du' 



(49) 



and ignoring terms 0(e), we can write 



I Aid 2 du/ = lim 



d 2 



uj'ujl iXe - 2 



0-ie 



Ah 2 c\uj 

iaj+xo+oo 



g 7ro;A x 



e -Vy-™\ d y 



lUl^XQ 



e -2/^^A d ~ da; / _ ( 5Q ) 



Let us now take the difference to the Xq — > limit 

/•oo /-co 72 

/ |/£,J a du/ - lim / |A^J 2 d^ = lim— °- 
Jo x °~*° Jo lZ A7rc i 



-e~ nuiX X 



uj'ujI iXe - 2 



0-ie 



ioj+xo+oo 



1UI+XQ 



e y dy 



\T(l-]u\)\' 



duj' . 
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Note that the expression inside the square bracket is well defined for all u/. We now 
replace the integration variable u/ by uj + , using its definition from above, and obtain 
(we omit e and e in the discussion of the following step) 



1 



-7tUj\ 







/■iaJ+xo+oo 








(/;< 













poo 




piw+xo+oo 








UJ / oj+ 2 




/ V 









-iwA -1/ 




e y dy 



The second term is finite (see Appendix) and, since we expect a divergent result, we 
can neglect it. The first term can be written as the integral from to oo minus the 
integral from to u. The latter integration is finite as well, since the difference of the 
T-functions vanishes as uo . — > 0. We are left with 



poo poo -1 

/ |/t J 2 da/ - lim / |/&J 2 da/ = finite + — « 



-e ~™ A x 



UJ 



-1 



iaj+xo+oo 



e y dy 



- |T(1 -iwA)| 5 



dcuj 



Now we neglect all finite terms and use the fact that the second term inside the square 
bracket is again (up to finite terms) the result of the xq — > limit (f45f ). We thus see 
that the divergent part of ([50]) can be written as 



OO -y 



e -™\ x 



x— j< piuj+XQ+oo p — iuj^XQ+oo 

iAe-1 / -\u)\ j / -y~-iw\ 



+ 



e 



dy / e-^ A dy d^ + . (51) 

'0— it Jiuj+xq J —ilj^xq 

For positive e we can replace the contour parallel to the positive real axis by an 
integration along the negative imaginary axis. Then we substitute e z = i— and have 

"OO 



r i 



-e~ nuiX X 



UJ 



fXQU]e z +oo 



e y dy 



xouie z * +oo 



(52) 



y iw V 5 dy dz . 

'ln(^) J xouje z J xoue z * 

We can now take the limit e — > and see that, since the integration is along the real 
axis, our integral can be written as 



poo poo 

/ |/td 2 du/ = Km / e-^ z h{z) dz 

J0 £ ^°J-oo 



with 



~iru)\ poo 



h{z) :-. 



4tt 2 lu 



y 



-\w\ -y 



e y dy 



y iujX e~ y dy , 



(53) 



(54) 



' xou)e z J XQLoe" 1 

for any complex z. For real z, as we integrate over on the r.h.s. of Eq. (|53|) , this yields 
the quantity we search for; h(z) is identical to the product of the inner integrations in 
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(|52|) for real z. This function is holomorphic in z. It has an upper bound B for any real 
z in [-R, R). Thus 



h(z)dz 



R 



< 2RB 



and, because the function is holomorphic, 



h(z)dz 



d 



< 2RB and 



h(z)dz 



< 2RB 



(55) 



(56) 



where the path 7^ is the half-circle of radius R in the lower complex plane and 7^ the 
half-circle of radius R in the upper complex plane. From this we find 



lim 




h(z) 


R—rOO 




z 2 




lim 


/ 


h(z) 


R^oo 




z 2 



dz\ 



dz\ 





. 



Our integral of interest can now be expressed as 

1 



h(z) dz 



-d z 
A 2 



iXez 



h{z) 



dz , 



(57) 
(58) 

(59) 



which allows us to make use of the residue theorem. Now for e > we may close the 
contour in the lower plane. To compute the integral we shift the pole into the lower 
plane. This yields 



Z 2 



With this we now get 

f°° 9tt 
/ |^,J 2 d^'= lim-h(0)5(e) 
Jo A 



lim - 



5(e) 



y 



-i^A -3/ 



e y dy 



y iujX e~My 



(60) 



(61) 



(62) 



' XQLU J XQUJ 

with the differentiation on the 5-function interpreted as 0.|] If we now put all pieces 
together, we finally find for the vacuum expectation value of the number density operator 
N(u) [in case (1)]: 

-tvujX 



l/t J 2 dJ 
N(uo) 



2tt\uj 



5(e)\T(l-iuj\,xouj)\ 2 + finite 
|T(1 — Xquj) | 2 . 



V 2tt\uj 1 

The spectral energy flux through a 2-dimensional surface p(u) then reads 



r(l + iu\, xquj) 



(63) 



(64) 



4tt 3 A 

We have made several assumptions, but some of them may be relaxed. First, we 
assumed that the velocity v of the mass shell is constant. In a realistic scenario we would 

S It is an even function. 
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expect a high velocity at the beginning, decreasing towards a constant value (inertial 
motion), in contrast to a collapse where the velocity would increase. Secondly, we 
approximated around the space-time points that cause the maximal amount of particle 
production. Therefore, our result is only correct in the vicinity of u = 0. To obtain a 
result for later times, we would have to approximate around other points. This gives 
rise to different values of the parameters A and Xq, which actually become functions of 
R and v(R). These functions then are [cf. fl§7])]: 

w , 2(1 + v)Md' 

\(R) = K ' and 

vd A 

x(R) := - A , with 

vd 



d'{R) 



d dV 



d-f(R) dv 



d(R) := ^ +d'{R)[l- 1 {R)} , (65) 

dV 7(i?) 

from which we obtain a spectral energy flux that depends on the actual radius and the 
actual velocity of the mass shell: 

p(u,R)= |r[l + i^A(/?),x(fi)]| 2 . (66) 

The quality of this generalization mainly depends on the assumption that the change of 
velocity is small, so that ^-contributions can be neglected. Should the acceleration of 
the mass shell get too high, the patches in which the velocity can be treated as nearly 
constant would get too small and effects from the boundaries would grow important. 



5. Discussion 



The particle production of an expanding mass shell depends on the velocity of the mass 
shell u, its mass M, and its radius R. In the spectral energy flux, these three quantities 
enter only through the functions A(i/, M, R) and x(u, M, R). To get a feeling for the 
behaviour of the function p(uS) we will first turn to a discussion of these quantities. In 
Table 1 we provide the values of A and x for some limits. The limit v — > 1 together with 
R — » 2M reproduces Hawking's result. 

In the Hawking limit we have A — > k, the so-called surface gravity. Additionally 
x — > 0, and we obtain exactly the expected Planck spectrum (cf. fl6"E|)). Figures |3] and f| 
show 1/ A as a function of R for three different values of R and x as a function of v for 
different values of R. We normalize all quantities to the Schwarzschild radius R$ = 2M; 
1/A decreases rapidly with decreasing velocity. As is shown below, this implies a decrease 
of the spectral energy flux and a shift of the maximum towards smaller frequencies as 
the velocity decreases; A is less sensitive to changes of the radius at high velocities. 
Particle production expected, at large radii, because x increases with large R 

and small velocity. 
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Figure 3. The function Rs/X(R) for different velocities v. 
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Figure 4. The function x(v)/Rs for different values of the radius R. 
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Limit 


A 


X 


v -> 1 


4M 


2(R- 2M) 


v -> 


oo 


oo 


i/ -> 1 A i2 -> 2M 


AM 





i? — >■ oo 


V 


oo 


M — > 





u 


i?^ 2M 


4M 






Table 1. Some limiting cases for the functions that enter the expression for the 
spectral energy flux through a 2-dimensional surface. 

5.1. The spectral energy flux 

We now study the behaviour of the spectral energy flux ( |S4]) as a function of X/Rs 
and x/Rs- Figures ^| and || show p(u) for different values of x/Rs, and for different 
values of X/Rs with the second parameter fixed at a representative value each time. 
Let us first turn to Fig. |5|. When x = we find a Planck spectrum with temperature 
T = ^TtksX)" 1 . For x/Rs > the spectrum is non-thermal. It is interesting to note 
that in the non-equilibrium situation more hard particles and less soft particles are 
produced at a fixed value of A. 

In Fig. ^]we see, as expected, that an increase in X/R$ causes a shift of the maximum 
towards smaller frequencies. Furthermore the whole curve flattens quickly for increasing 
X/Rs- 

Remember that, according to the considerations presented in the last section, an 
increasing radius (which changes A and x) dominates the evolution of the spectral energy 
flux during the expansion. This can be seen in Fig. [7|. 

5.2. The energy flux 

Finally we come to the emitted energy flux through a 2-dimensional surface, e = 

In Fig. [8] we see the energy flux as a function of radius for different values of v. At 
large values of R/Rs no particles are produced, since the mass shell moves inertially 
(the tidal forces go to zero). However, more interesting is the maximum of the curve 
that appears for ultra-relativistic velocities v. This reflects the fact that two effects are 
working against each other during the expansion. The probability of particle production 
depends on the strength of the tidal forces acting on the vacuum locally. The tidal 
forces decrease monotonically as the mass shell expands. Once particles have been 
created, they have to climb out of the gravitational potential well of the mass shell itself 
(gravitational redshift). This effect reduces the energy of the produced particles. At the 
Schwarzschild radius the energy would be reduced to zero. For an ultrarelativistic mass 
shell the tidal forces change very quickly and thus the particle production is large. Close 
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Figure 5. The spectral energy flux through a 2-dimensional surface p(u>) for different 
values of x/Rg at fixed X/Rs = 2.05. 




Figure 6. The spectral energy flux through a 2-dimcnsional surface p(u>) for different 
values of X/Rs at fixed x/Rs = 1. 
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Figure 7. The spectral energy flux through a 2-dimensional surface p{to) at different 
moments of time, respresented by three positions R/Rs, for the constant velocity 
v = 0.95. 




Figure 8. The energy flux through a 2-dimensional surface for different values of the 
velocity v. 
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to the Schwarzschild radius however, this is counteracted by the gravitational redshift, 
giving rise to a maximum slightly away from i?s- 



6. Summary 

We computed the spectrum of particles produced by the time-dependent gravitational 
field of an expanding mass shell. This system is — as long as no horizon forms — invariant 
under time-reversal, so that we can describe an expanding shell as well as a collapsing 
shell before the formation of a horizon. If the horizon is formed, the system gets quasi- 
static for the asymptotic observer and the effect is "frozen" . In this way we can reproduce 
the result of Hawking radiation. 

In the situation of an expanding mass shell the spectrum is non-thermal. Our result 
might be of relevance for the discussion of evaporating black holes. In the situation of 
a collapsing mass shell, we obtain the deviations from the Planck spectrum of Hawking 
radiation, which is produced just before the horizon forms. 
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Appendix A. 



To show that 



UJ 



io;+xo+oo 



-Ud+ XO + OO 

y-^ e -Vdy-\r(l + h}\)\ 2 



do;+(A.l) 



(A.2) 



y iu}X e- y dy 

is finite (for uj > 0), it is sufficient to show that 

y iujX e~ y dy < const., 

which means that it has an upper bound for uo + — > oo since the second term in ( |A.1| ) 
does not depend on u + anyhow. We are then left with an integral over uj+ 2 that is 
finite. If we succeed in showing that (|A.2| ) is true, we know that the first term does not 
increase faster for large uj + and therefore is finite too. It is 



lU+XQ+OO 



y iujX e~ y dy 



IX 



iuj+xa+oo 



3 iwAln|j/|-£jAaxg(y) -y 



< sup 

j/67 



i^Aln|j/|-wAarg(j/) 



e y dy 

iu;_|_oo 



e- y dy 



XQCO+ 



(A.3) 



(A.4) 



7 being the path of integration in the upper right quadrant of the complex plane and 
therefore is arg(y) e [~7r, 2ti\ on 7. This yields 



y iujX e- y dy 



1XQUJ+ 



< e 



1XQUJ+ 



e y dy 



(A.5) 
(A.6) 
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Putting this together we have 




y dy <e 



(A.7) 



which shows that (|A.1|) is finite, as claimed. 
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